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Abstract. We compare the current effects of rotation in stellar evolution to those of the magnetic field created by 
■ the Tayler instability. In stellar regions, where magnetic field can be generated by the dynamo due to differential 

' rotation fSpruit l!jU02ll . we find that the growth rate of the magnetic instability is much faster than for the thermal 

instability. Thus, meridional circulation is negfigible with respect to the magnetic fields, both for the transport of 
angular momentum and of chemical elements. Also, the horizontal coupling by the magnetic field, which reaches 
^ ' values of a few 10^ G, is much more important than the effects of the horizontal turbulence. The field, however, is 

not sufficient to distort the shape of the equipotentials. We impose the condition that the energy of the magnetic 
field created by the Tayler-Spruit dynamo cannot be larger than the energy excess present in the differential 
rotation. This leads to a criterion for the existence of the magnetic field in stellar interiors. 

Numerical tests are made in a rotating star model of 15 Mq rotating with an initial velocity of 300 km-s^^. We 
' find that the coefficients of diffusion for the transport of angular momentum by the magnetic field are several 

, orders of magnitude larger than the transport coefficients for meridional circulation and shear mixing. The same 

applies to the diffusion coefficients for the chemical elements, however very close to the core, the strong ^-gradient 
reduces the mixing by the magnetic instability to values not too different from the case without magnetic field. 
We also find that magnetic instability is present throughout the radiative envelope, with the exception of the very 
Q ' outer layers, where differential rotation is insufficient to build the field, a fact consistent with the lack of evidence 

)-i I of strong fields at the surface of massive stars. 
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1. Introduction are considering here following Spruit fl9y9' "20021. Sect. 

3 compares the characteristic times of meridional circula- 
■ - ' The inclusion of new physical effects in stellar evolu- ^-^^ magnetic field instabilities. Sect. 4 considers 
tion greatly improves the comparisons with observations. ^^^^ happens to the horizontal turbulence in presence of 
About twenty years ago the impact of mass loss by stellar j^agnetic fields. Sect. 5 shows a new physical limit on the 
wmds on the evolution was found to be large. However, occurrence of the magnetic field in rotating stars. Sect. 6 
some significant discrepancies remained and the inclusion gj^^g ^^^^ numerical values on the size of magnetic and ro- 
of rotation has enabled substantial progress in the com- tational effects in the case of a 15 Mq star. Sect. 7 presents 
parisons with observed chemical abundances, with number conclusions 
counts and with chemical evolution of galaxies (cf. Langer 
et al. ll999l Maeder & Mevnet .2001..) . Magnetic field is com- 
ing next, but certainly not the last, in this series of effects 



2. Basic properties of the magnetic field 



which may influence all the outputs of stellar evolution. Let us collect here some basic expressions and concepts 
In this work, we focus mainly on the relative impor- we need below. Spruit H1999L 1^028 has shown that the 
tance of the effects of the magnetic field and of rotational magnetic field can be created in radiative layers of stars 
instabilities to try to determine which effects can be let in differential rotation. Even a small toroidal field is sub- 
aside and what must be considered as a priority. Sect, ject to an instability (called Tayler instability by Spruit), 
2 summarizes the main effects of the magnetic field we which creates a vertical field component. Differential rota- 
tion winds up this vertical component, so that many new 
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lines become progressively closer and denser in a star in a 
state of differential rotation and therefore a much stronger 
horizontal field is built. This is the dynamo processs de- 
scribed by Spruit. Tayler instability is a pinch-type in- 
stability. As shown by Spruit H1999|l . it has a very low 
threshold and is characterized by a short timescale. Also 
it is the first instability to occur. The magnetic shear in- 
stability may be present, but it is of much less importance 
(Spruit 1999). 

The instability occurs in radiative zones and two cases 
are considered by Spruit H1999I I2()()2(l depending on the 
thermal and /i~gradients through the associated oscilla- 
tion frequencies, 



m = 



and 



9l_ 
Hp 



(V - Vad) 



9^ , 
Hp 



(1) 



(2) 



The thermodynamic coefficients 5 and ip are defined as fol- 

lows S = - (l^)p^ and ^ = (flf^)^^p. The quantity 

Hp is the local pressure scale height. Following Spruit, we 
call case the case where the /i-gradient dominates over 
the thermal gradient, i.e. when Nf^ > N-r. Case 1 applies 
when the thermal gradient is the main restoring force, i.e. 
when < Nt. Let us point out that cases and 1 are the 
two simpler limits of a more general case, where a proper 
account of both thermal and /i effects would be made. We 
may wonder whether we miss some important physical 
situations with this simplification. Indeed, numerical tests 
show that in practice the intermediate situation between 
cases and 1 introduced by Spruit concerns a small but 
significant part of the star as shown in Fig.|21 In this part, 
cases and 1 give diffusion coefficients which have the 
same order of magnitude, therefore the general effect is 
correctly described. However, there is little doubt that in 
future the general non-adiabatic case has to be examined 
in detail. 

The growth rate of the magnetic instability is 



n 



(3) 



where to a is the Alfveri frequency, i.e. the frequency of 
magnetic waves. As shown by Pitts & Tayler H1985fl . the 



reduction factor ^ in Eq. (O results from the Coriolis 
force in a star rotating with angular velocity il. The Alfven 
frequency is 



B 



LOA 



(4) 



The magnetic instability works only if the unstable dis- 
placements do not loose too much energy against the sta- 
ble stratification. For this to be the case, the radial dis- 
placements (against the buoyancy force) must be small 
compared to the horizontal displacements. Taking r as a 



maximum for the horizontal displacements, this sets an 
upper limit on the radial length scale ZrO,i of the displace- 
ments, 



and 



r UJA 
N,, 



K 



(5) 



(6) 



where the indices or 1 refer to the two cases considered. 
The term K is the thermal diffusivity K = ^0^, where 
the other quantities have their usual meaning in stellar 
structure. There is also a minimal extent of the magnetic 
instability, below which it is quickly dissipated by mag- 
netic diffusivity. 



It 
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(7) 



where rj is the diffusivity of the magnetic field. The com- 
bination of Eq. {7)) with Eqs. Q and © leads to a min- 
imum value of the Alfven frequency for the occurrence of 
the magnetic instability in the two cases considered, 



(^) > 
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At 
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(9) 



These are the conditions in order that the magnetic field 
overcomes the restoring force of buoyancy. In addition in 
the second case, the effects of the thermal diffusivity de- 
scribed by K are also accounted for. The corresponding 
maximum radial dimensions of the instabilities are given 
by the above Eqs. (OEl- 

Spruit H2()()2I Eqs. 18, 19) considers the amplification 
timescale necessary to double the component B^ start- 
ing from the radial component B^ over the largest char- 
acteristic lengths defined by Eqs. (jsUHJ. He assumes the 
equality of the amplification timescale with the timescale 
for the damping by magnetic diffusivity over the above 
lenghtscales. In this way, he obtains the expressions for 
the Alfven frequency. In the first case, where A^^ domi- 
nates, this is 



(10) 



where g = Thus, we see that the Alfven frequency 

(a measure of the field amplitude) depends on the differ- 
ential rotation parameter and on the ratio of the angular 
velocity to the Brunt -Vaisala frequency. When thermal 
diffusion is accounted for (A'^ negligible), one has 



At 



K 



r2 A^T 



(11) 



there the thermal diffusivity also intervenes. A few other 
recalls will be made when necessary. 
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3. Growth rate of the magnetic instability with 
respect to meridional circulation 

A major question arises concerning a rotating star with a 
magnetic field. What happens to the meridional circula- 
tion in presence of the magnetic field of the Tayler-Spruit 
dynamo? Basically, meridional circulation occurs because 
thermal equilibrium cannot be achieved on an equipoten- 
tial inside a rotating star. Thus, we may wonder whether 
the horizontal breakdown of thermal equilibrium in a ro- 
tating star can be compensated by a magnetic stress on 
an equipotential. 

Let us define a velocity C/magn characterizing the ra- 
dial growth of the magnetic instability. In the two cases 
and 1, we consider the ratio of the appropriate maximum 
lengths given by Eqs. lOEJ to the characteristic time fi^^. 



Kiiagn, 0,1 — 0,1 f • 

For the case with iV^ > iVr, one gets 



^magn, 



(12) 



(13) 



Then, using the above expression (Eq. Il()|l for the Alfvcn 
frequency, we obtain 



t^magn, = f?" r ( 



(14) 



The magnetic instability grows very fastly with the angu- 
lar velocity and the differential rotation parameter q. In 
case 1 with iV^ < iVx, we get from Eq. ((T^ with Eqs. © 
and © 



^inaen, 1 — ^ 



(15) 



Using Eq. 1)11(1 for the Alfven frequency, one has for the 
growth velocity of the magnetic instability 



q r Q 



K 



(16) 



This velocity also grows with Q and q, but less than in 
the case 0, because the radiative diffusivity reduces the 
dependence of the Alfven frequency with respect to rota- 
tion parameters Q and q. 

These velocities have to be compared with the radial 
component of the radial part of the velocity of the merid- 
ional circulation U{r), as given by Maeder & Zahn ( 1 998,. 
Eq. 4.38). If the circulation velocity would largely domi- 
nates, this would mean that the magnetic field has effects 
which are too weak to infiuence the meridional circulation 
and the circulation would develop as usually supposed. If 
on the contrary, one has ?7magn, o,i >> U(r), this means 
that the magnetic instability develops much faster than 
the thermal instability at the origin of meridional circula- 
tion. If so, this means that the thermal instability created 
by rotation on an equipotential will interact firstly with 



the magnetic field. Tayler instability, which has the short- 
est timescale, may possibly develop and create a magnetic 
field which will introduce some stress horizontally on the 
equipotential. Detailed calculations must be done with ac- 
counting for the effects of the magnetic field in the equa- 
tion for the entropy conservation at the basis for the cal- 
culations of meridional circulation, (this has been made 
for the effects of horizontal turbulence on the meridional 
circulation by Maeder & Zahn|lg985. 

For now, we assume in this case, in the whole region 
where C/magn, o,i >> U{r), that the usual circulation ve- 
locity must be set to zero. This working hypothesis is 
largely confirmed by the comparison below (cf. Tables 1 
and 2) of the velocities U (r) of meridional circulation to 
the velocities ?7magn,o,i characterizing the growth of the 
magnetic instability, the last ones being 4 to 7 orders of a 
magnitude larger than the first one. 

Before looking more to the numerical values, we note 
that some caution has to be taken so that the compar- 
isons of C/magn, 0,1 ^ud C/(r) are done in a consistent way. 
The expression for the transfer of angular momentum by 
circulation and diffusion D is (see Zahn fTO^ Maeder & 
ZahnHMl 



1 d 



J.2 Qj, [P Qj, 



(17) 



where Q{r) is the average angular velocity on the equipo- 
tential. From the second member of this expression, we 
see that to any diffusion coefficient D, one can define an 
associated velocity Ub 



D dn D 
= 5 - — = h-q 
It or r 



(18) 



The factor 5 results from the integration of the angu- 
lar momentum over the star, which leads to the above 
Eq. (|17|) . In order to be correct, the comparison of ve- 
locities must be made over the corresponding quantities 
as they appear in Eq. H17|) for the angular momentum 
transfer, which in turn implies Eq. I|18|) . Now, we take 
the coefficients of magnetic diffusion rjQ and rji , which are 
obtained by replacing the inequalities of Eqs. (|H1 © by 
equalities. This means (cf. Spruit ,2002J that we consider 
the cases of marginal stability. Now, we can associate in- 
tegrated velocities to these diffusion coefficients thanks to 
Eq. H18|) and we find that they are just equal to the ve- 
locities C/magn, and C/magn, 1 givcu abovc multiplied by a 
factor of 5. These velocities characteristic of the magnetic 
instability given by Eq. (|18|l are those to be compared to 
the velocity of meridional circulation, as it appears in the 
equation for the transport of angular momentum. As we 
shall see below, even if we omit these various considera- 
tions about Eq. (|17|l . which result in the above factor of 
5, the conclusions on the order of magnitude will be the 
same. 
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Below in Sect. 6, a star model with an initial mass of 15 
M0, a composition given by X = 0.705 and Z — 0.02 and 
an initial rotation velocity of 300 km-s~^ has been com- 
puted. The prescriptions for rotation are those in Maeder 
& Meynet H2000|l . Tables 1 and 2 show the main param- 
eters when the central H-content is = 0.60 at an age 
of 4.17 • 10^ yr. We see that the velocities 5 ?7magn, 0,1 are 
3-10^ to 3-10^ times larger that the velocity Uir) of merid- 
ional circulation. This shows that, even if present, the ve- 
locity of meridional circulation is negligible with respect 
to the corresponding velocity characterizing the transport 
of angular momentum by the magnetic field. Thus the 
question whether or not meridional circulation appears in 
the presence of magnetic field is of little relevance, since 
meridional circulation would anyway be totally negligible 
in comparison of the effect of magnetic instability. The 
same remark arises when one compares the magnetic dif- 
fusion coefficients for the angular momentum as given in 
Eqs. (@ni 1^ with the corresponding expression 
for meridional circulation. 

Thus, we suggest that in the stellar regions where mag- 
netic field is present, and only there (see Sect. 5), the 
meridional circulation may be neglected with respect to 
magnetic field effects for the transport of angular momen- 
tum. As is shown in Sect. 4 below, the above conclusion 
also applies to the transport of the chemical elements. 

4. The importance of magnetic stress vs. 
horizontal turbulence 

4.1. Horizontal coupling 

The turbulence in rotating stars is highly anisotropic and 
has a strong horizontal component, described by a diffu- 
sion coefficient Dh, because vertically the thermal gradient 
stabilizes turbulence. This horizontal turbulence strongly 
reduces the horizontal differential rotation, so that rota- 
tion varies only radially. The rotation is therefore said to 
be "shellular" (Zahn 119921 f2 uniform at the surface of 
isobaric shells) . The coeflicient Dh also plays a role in the 
mixing of the chemical elements and meridional circula- 
tion (Maeder & Zahn HM^ . A first expression for Dh was 
given by Zahn H1992|l . Recently another expression of this 
coefficient has been obtained (Maeder | 



Dh^Ar [rn{r) V {2V - aU) 



(19) 



with A « 0.10, V{r) the horizontal component of merid- 
ional circulation and a = ^ '^dinr^ ■ Typically, this new 
estimate of Dh is of the order of 10^"'^ to 10^^ cm^-s^^in 
a massive star, i.e. typically 10^ to lO'^ times larger than 
the coefficient previously estimated. Recent developments 
by Zahn (private communication) and entered into model 
calculations by Palacios (private communication) confirm 
this higher values of Dh. 

The question is now: what happens to this horizontal 
turbulence, even if it is much larger than initially sup- 
posed, in the presence of the magnetic field? According 



to Spruit H2002fl . the Tayler instability and the associated 
dynamo leads to horizontal field components in cases 
and 1, 



B^,o = (4 7rp)2rg 



B^,i = (4 TT p)2r 17 



n \^ / K 



(20) 



(21) 



For the numerical model of Tables 1 and 2, we find an 
horizontal field 



= 4.35 • 10^ G at M,/Mq = 6.57 , 
= 2.80 • 10^ G at M,./Mq = 11.04 



(22) 
(23) 



These fields are high with respect to our current stan- 
dards, but the magnetic pressure. 



Pr 



magn 



(24) 



is small with respect to the total pressure. One has typ- 
ically a ratio Pmagn/P = 1-5 • 10~6, 2.7 • 10"^ ^he 
above examples. Thus, the magnetic field generated by 
the Tayler-Spruit mechanism does not affect the shape of 
the equipotentials. The ratios B,-/B^ of the radial to the 
horizontal field components are given by Eqs. (21) and 
(23) by Spruit H2002|l and the values are of the order of 
10-3 to a few 10-2. 

Now, the question is how does the horizontal coupling 
due to the magnetic field i?ip,o,i compares with the hori- 
zontal turbulence characterized by Dh. Spruit l|2002l Eqs. 
31, 32) has given an estimate of the coefficient of the verti- 
cal coupling by magnetic field. However, we cannot use it 
here, because the field is very anisotropic. The coefficient 
Db,, for the horizontal coupling must be much larger than 
the vertical one. At low rotation, we suggest to take the 
following estimate. 



2 ^ B^ 



(47r p): 



(25) 



At high rotation, the growth rate is also reduced by the 
Coriolis factor u!\/^ (as suggested by Spruit in a private 
communication), thus one has then 



Bt 



(26) 



47r p f2 

Interestingly enough, this coefficient is of the order of 10^'* 
cm^-s-^ in the above examples. Thus, the horizontal mag- 
netic coupling is about 10'^ larger than the horizontal cou- 
pling by turbulence, even when we take the value given 
by the new coefficient Dh in Ea. p9|) above. Therefore, 
our conclusion is that in the presence of Tayler-Spruit 
dynamo, we can anyway neglect the coupling due to hor- 
izontal turbulence with respect to the coupling insured 
by the horizontal component of the magnetic field. More 
likely, we suspect that the horizontal turbulence may not 
develop in this case or occur in a very different and limited 
way. 
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4.2. Remarks on further consequences: the case of Dcs 

In models with rotation and without magnetic fields, the 
combined effect of meridional circulation and horizontal 
turbulence may be treated as a diffusion with a coefficient 
usually called Z?off (cf. Chaboyer & Zahn lTM^ . 



D, 



off 



r U{r)Y 
30 



Now, since I?h is negligible, we may wonder what hap- 
pens to the transport of chemical elements by meridional 
circulation. As well known, it is in general not correct to 
express an advection as a diffusion, but let us do it here 
just for the purpose of comparing orders of magnitude. 
The diffusion coefficient we may associate to the merid- 
ional circulation is in this case. 



\rU{r)\ 



(28) 



This must be compared to the coefficients Dchem, o, i, 
which are obtained by imposing equality in Eqs. (|H1 El 
and taking -Dchom, o, i — Vo,!- The equality means that 
one considers the case of marginal stability as determin- 
ing for the transport. Then, in the obtained relations, the 
Alfven frequency is expressed with the help of Eqs. 1)101 
[TT|) . In this way. Spruit obtains l|2nn^ Eqs. 42, 43), 



D. 



chcm, 



and 



chcm, 1 



"'111 



K 



N'T 



(29) 



(30) 



In the numerical example of Tables 1 and 2, we find ratios 
= 0.11, 3.2 • 10"^, 1.9 • lO"'' at the 3 layers consid- 
ered, (taking in each case the appropriate expressions for 
the case and 1). Thus, we see that through most of the 
star the circulation (if it exists!) would have a negligible 
effect for the transport of chemical elements. Only close to 
the core, the circulation might represent about 10 % of the 
transport by the magnetic instability. This is small any- 
way and owing to the profound doubts expressed above 
on the occurrence of any circulation, we consider that any 
transport of chemical elements by meridional circulation 
can generally be neglected in the presence of Tayler-Spruit 
dynamo. 

5. An energy condition for the magnetic field 

5.1. Energy condition 

In a radiative zone, the magnetic field created by the 
Tayler-Spruit instability arises from differential rotation. 
Therefore, we must impose the condition that the en- 
ergy in the magnetic field cannot be larger than the ex- 
cess energy in differential rotation. Tayler 1)1973(1 and Pitts 
and Tayler fl986) have studied the conditions for the ap- 
pearance of the magnetic instabilities, which have very 



short growth times. The above condition is different, it 
is an energy condition, which must be satisfied on longer 
timescales. Due to the magnetic diffusivity, the field once 
created tends to disappear. Taking the values of Dchcm = V 
given for example in Fig. 6 (see also Table 1 and 2), which 
are between 10^°'^ and 10^^ cm^-s~^ over most of the star 
except very close to the convective core, we find that the 



(27) timescale 



for the diffusion of the magnetic field 



is of the order of 2 -lO"^ to 8 TO^ yr., which is short with 
respect to the stellar lifetimes. This means that the field 
created by Tayler-Spruit dynamo will exist only if the en- 
ergy of the magnetic field is continuously replenished by 
differential rotation during the MS evolution. Therefore, 
we need to apply the above condition. 

For a magnetic instability with a displacement of am- 
plitude ^, the kinetic energy Eq by unit of mass is 



(31) 



The excess energy Eq in the differential rotation is the 
difference of energy between the existing flow with differ- 
ential rotation and a flow with an average rotation over the 
considered radial distance r. Let us consider a horizontal 
velocity field W{r). Over a vertical distance dr, the energy 
excess dE^ over the extent of the magnetic instability is 



1 



1 



dEn = - [W^ + iW + dW)^]---2iW 



dW 



dW^ 



.(32) 



If and /h are the radial and horizontal components of 
the displacement ^ of the magnetic oscillation (cf. Spruit 
12002V the energy excess over the displacement ^ can be 
written. 



E. 



1 



fdWV 

4 yd^j 



(33) 



The horizontal velocity W{r) — rsin-i? fi, where d is the 
colatitude. Thus, one has 



dW 
dr 



dn d\nn . 

r sm V —r- — s I sm v — = —il sm v q 



dr 



d\nr 



The condition Eq > E^ leads to 
A /h 



(34) 



(35) 



This would apply to a given colatitude "d. This equation 
as it stands means that the reservoir of available rotational 
energy is larger at the equator, while we know (cf. Spruit 
I1999|l on the other side that Tayler instability is stronger 
away from equator. Now, for the ratio (j^) of the vertical 
to the horizontal displacement, we could take for example 
in the outer layers, where differential rotation is generally 
small, the value given by Eq. © above (case 1). This 
promptly leads to the following criterion with account of 
Eq. inj, 



l9l>3 ^ 



r^Ni 
K 



(36) 
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This equation ignores the geometry of the field. However, 
we have to consider carefully the geometry of the problem 
in 2 specific respects: 

1.- We must account for the fact that in the model of shel- 
lular rotation the energy of rotation is not a local quantity 
depending on (r, z?) , but on r only. The physical reason in 
usual rotating models is the strong horizontal turbulence 
(cf. Zahn llOi?^ . In the present models, the horizontal mag- 
netic coupling is even stronger as seen above in Sect. 4.1, 
so that shellular rotation is a valid assumption here. In 
such a case, the average stellar structure of the rotating 
star very well corresponds to the structure at a colatitude 
iD given by the root of the second Legendre polynomial 
P2{cos'd) = 0, (this has been verified in a recent work, 
Maeder "SOOl ). Thus it is appropriate to consider for the 
differential rotation in Ea. H34|) on a given equipoten- 

tial the average value (^)^ = | ^^q^. 
-2. The geometry of the field is also particular (cf. Spruit 
119991 12()()2|l . It consists in stacks of magnetic loops con- 
centric with the rotation axis. The main component of the 
displacement due to Tayler instability is perpendicular to 
the rotation axis. This means that at colatitude 1}, the ra- 
tio (j^) behaves essentially as tanz?. Since the polar caps 
are most unstable, while the equatorial regions are not, 
it is thus clear that in polar regions one has a ratio (j^) 
smaller than 1. However, in 1-D models as here we must 
consider the significant average for the whole range of co- 
latitudes. The field behaves as B^p sa sin-iJcosi? (cf. Spruit 
119991 Eq.35) and the Tayler instability develops only for 
< 7r/4, therefore it is necessary on a given isobar to 
consider colatitudes 1) smaller or at most equal to tt/A. If 
we take this last value as the limit, this gives the upper 
bound tanz? = (f-) < 1. 

From these two geometrical remarks, we obtain the 
necessary condition for the existence of a magnetic field 
generated by the Tayler-Spruit dynamo as follows, 



|g| 



(37) 



It means that the degree of differential rotation q must at 
least be larger than a/3 times the ratio ^ of the Alfven to 
the rotation frequency, in order that there is enough en- 
ergy in the differential rotation to allow the Tayler-Spruit 
dynamo to operate and build a magnetic field. We insist 
that this is only a necessary condition. If this condition is 
not realized, there is certainly no magnetic field created by 
the dynamo. Further works may perhaps lead to an even 
more constraining condition. The numerical factor, here 
\/3, may depend on the exact geometry of the magnetic 
displacements in a rotating star. We note also that if there 
are several types of instabilities generated by differential 
rotation, the available energy given by Eq. H33(l would in 
some way be shared between the instabilities. However, as 
mentioned above, the Tayler instability is the main one 
and shear instabilities appear negligible in comparison. 

We can go a step further, since the Alfven frequency 
LOA is a function of rotation and differential parameter 



\q\. For the case 0, where the /^-gradient dominates, the 
ratio ^ is given by the above Eq. ifTUI) . Thus, the above 
condition H37|) becomes in this case 



n 1 

— < ^ = 0.5774 

V3 



(38) 



This is the necessary condition in order that a magnetic 
field may develop from differential rotation in case 0. At 
first glance, this condition may look strange, since it means 
that f2 must be smaller that some value. The reason is that 
to A grows like q , while the upper limiting oja expressed 

by Eq. ^ goes hke ^ n. Thus, if n would be too big, 
the actual to a would overcome the critical value. In the 
numerical examples of Tables 1 and 2, case is relevant at 
the edge of the core at M,/Mq = 6.57. There is 0.11, 

(typically lies between 0.1 and 0.15 throughout the 
star). This is smaller than 0.5774 and thus the magnetic 
field can be present in these layers. 

We now consider case 1 with thermal diffusion, the 
Alfven frequency is given by Eq. (|ll|l . With the above 
condition 13 7|) . one obtains 



> 



K 



(39) 



This means that the differential rotation parameter \q\ has 
to be large enough to be able to generate the magnetic 
field. In the numerical examples in Sect. 6, the condition 
is not satisfied in the very outerlayers, which have a too 
weak differential rotation. 

5.2. Collection of formulae and recipes 

Let us collect here the various expressions we have for the 
diffusion coefficients by the Tayler-Spruit dynamo in ra- 
diative zones. Case applies when N^^ > N^. Magnetic 
field is present only when the criterion given by Eq. H38|l 
is satisfied. Then the diffusion coefficients for the trans- 
port of the angular momentum and chemical elements are 
respectively. 



Da.ngO = r^^q^ ( 



Ddacma = r'^^q^ ( 



(40) 



(41) 



Case 1 applies when < Nt, then account is given to 
the thermal diffusivity. Magnetic field is present only when 
the criterion given by Eq. (|39|1 is satisfied. The diffusion 
coefficients for the angular momentum and chemical ele- 
ments are respectively 



-Dchcmi = r-fl \q\ ( — 



1/2 



n 



K 



3/4 



1/2 



K 



3/4 



(42) 



(43) 
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Fig. 1. Internal H-profile in the 15 Mq test model with 
Willi = 300 km/s. The model is at an age of 4.0-10^ yr. This 
is the reference model in which we examine the properties 
of the magnetic field in detail. 




2x10 



r/Re 



Fig. 2. The oscillation frequencies in the model of Fig. ^ 
N'^ is indicated by a continuous line and is given by a 
dashed line. 



Normally, these coefficients are larger than those which 
would be obtained in case with Nt instead of N^, be- 
cause the account for thermal effects reduces the buoyancy 
force which opposes to the magnetic instability. However, 
as noted by Spruit (,,2002|i) . it may happen in some cases 



that these coefficients with index "1" are smaller. As noted 
by Spruit H2002|l . this is an artefact from the simplifica- 
tion introduced by considering only the 2 limiting cases 
and 1. Spruit suggests to introduce an interpolation for- 
mula depending on q. We hesitate to do so, because this 
extra-dependence on the differential rotation parameter q 
is unphysical, since the interpolation should rather depend 
on the thermal and magnetic diffusivities K and rj. Thus, 
the suggested treatment might introduce spurious effects 
in some evolutionary stages. The general case where ther- 
mal effects and /i~gradient are accounted for needs to be 
worked out in future. For now, we prefer to do the follow- 
ing by considering the coefficients 



irlP 



D. 



chcmlP 



r^nq" 



]Vt 



(44) 



(45) 



These are the same equations as in case 0, but with iVr 
instead of Nfj_. This means that we are considering only 
the restoring force of the thermal gradient and that we 
are ignoring the non-adiabatic radiative losses. 

For case 1, we compare the coefficient with indices "1" 
and "IP" and must take the larger ones. From Fig. lO 
below, we see that the coefficients "IP" may be larger 
than the coefficients "1" in some parts of the star. Of 
course, we have also to test whether the magnetic field 
can be created from differential rotation. For such zones 
of case IP, the criterion for the existence of the magnetic 
field is evidently the following one 



O 1 

— < ^ = 0.5774 , 



(46) 



and the appropriate test has to be made in the concerned 
layers. 

6. Numerical applications to stellar models 

6.1. Magnetic vs. rotational transports of angular 
momentum and of chemical elements 

We calculate evolutionary models of a 15 M0 star with 
initial composition X — 0.705 and Z — 0.02. The physics 
of the models (opacities, nuclear reactions, mass loss rates, 
treatment of rotation, increase of the mass loss rates with 
rotation, etc..) are the same as in Meynet & Maeder 
(2003). We compute the MS evolution of a model with an 
initial velocity of 300 km-s~^, which leads to an average 
velocity during the MS phase of about 220 km-s~^, which 
corresponds to the observed average rotation velocity. We 
now consider in detail the properties of a particular model 
with rotation to see the various coefficients and criteria 
characterizing the growth of the magnetic field. We take 
the model at an age 4.17 • 10® yr with a central H-content 
Xc — 0.571. The H-profile inside the star is illustrated 
in Fig. n and the oscillation frequencies and in 
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Fig. 3. The distribution of the angular velocity in the ref- 
erence model of Fig.l with rotation (continuous line), 




I 

I 
I 



I 1 ' '— 1 ' ' ' 1 ' ' ' 1 

2 4 6 

Fig. 4. The diffusion coefficients Uango (continuous line) 
and -Dchcmo (dashed line)corresponding to case 0. As dis- 
cussed in the text, these coefficients apply only in the ris- 
ing part between r/R© = 1.32 and 1.59. 

Fig. [3 We see that case applies between t/Rq — 1.32 
and 1.59, which corresponds to mass coordinates 5.53 and 
7.70 M0 respectively. The internal profile ^{r) is illus- 
trated in Fig. El 

The diffusion coefficients Dango and I?chcmO are illus- 
trated in Fig. ^ These coefficients grow very fastly at the 



Fig. 5. The diffusion coefficients for angular momentum 
i'angi (continuous line) and -Dangip (dotted line) . The dif- 
fusion coefficients for chemical elements I?chcmi (dashed 
line) and -Dchcmip (long dashed line). In each case, the 
largest diffusion coefficient has to be taken. 



edge of the core since decreases very fastly there and 
the dependence in the ratio goes like the power 4 

and 6 for the two coefficients respectively. These coeffi- 
cients apply only in the rising part between t/Kq = 1.32 
and 1.59 as indicated. Above this value, the main restoring 
force is no longer the /^-gradient, but the stable tempera- 
ture gradient. 

In case 1, where thermal gradients dominate, the dif- 
fusion coefficients are illustrated in Fig. [51 We see that 
for the transports of angular momentum and of chemi- 
cal elements, the coefficients with indices "1" dominate 
in the external parts of the star above t/Rq — 2.69. 
However, we also notice that in a sizeable region, i.e. be- 
tween t/Rq = 1.59 and 2.69, the coefficients with indices 
"IP" dominate over those with "1". This occurs, as ex- 
pected, at some limited distance above the edge of the con- 
vective core, at the place where the /i-gradient becomes 
small enough, but is still different from zero. In these re- 
gions, more general developments having case and 1 as 
limiting cases would be a progress. We see that the dif- 
ferences between the two cases "IP" and "1" amounts to 
a maximum of 0.4 dex and 0.6 dex for the transport of 
angular momentum and chemical elements respectively. 
This is limited, but non negligible, and it may justify a 
further study of the physics of the general case. However, 
we notice that this difference is small when compared to 
the differences resulting from the inclusion of the magnetic 
field or not, which as shown below amounts to several or- 
ders of magnitude. Thus, we conclude that the present 
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Fig. 6. The figure shows the complete description of the 
diffusion coefhcient I?ang for angular momentum (contin- 
uous line) made according to prescriptions of Sect. (5.2). 
The complete description of the diffusion coefficient -Dchcm 
for chemical elements is given by the dashed line. The dot- 
ted line shows the shear diffusion coefficient -Dghcar from 
Maeder H1997|l . The thermal diffusivity coefficient K is 
represented by a long-dashed line. 

coefficients of diffusion need to be further improved, but 
they nevertheless describe correctly the main results of the 
inclusion of the magnetic field. 

Tables Hand 121 provides some useful structural param- 
eters, the diffusion coefficients and the velocities of merid- 
ional circulation and of magnetic instability at 3 locations 
in the reference model of 15 Mq with an initial velocity of 
300 km s~^ at an age 4.17 • 10^ yr. The three levels con- 
sidered illustrate the case 0, IP and 1 respectively. These 
Tables permit further quantitative analysis of the various 
terms intervening in the equations. 

Fig. |H| shows the comparison of the diffusion coeffi- 
cients due to the magnetic field compared to the diffusion 
coefficient Dshcar by shear instability in the rotating star 
and to the thermal diffusivity K . The transport of angular 
momentum by the magnetic field is 6-7 orders of magni- 
tude stronger than by shear instability in a rotating star. 
Similarly, as discussed in Sect. (4.2), the magnetic trans- 
port of angular momentum is also 4-7 orders of magnitude 
larger than by meridional circulation. Therefore, we con- 
clude that transport of angular momentum by magnetic 
field is totally dominating, if magnetic field is present. 

For the transport of chemical elements, the difference 
between the two diffusion coefficients amount to 3-5 or- 
ders of magnitude in favour of the transport by magnetic 
field. The difference is especially large in deep regions at 
some distance of the convective core. At the very edge 



Fig. 7. The difference — ^ when case applies, i.e. 
in the region between t/Rq = 1.32 and 1.59. According 
to criterion H38|) . if this difference is negative, magnetic 
field is created. Thus, we see that in the region where the 
//-gradient dominates, magnetic field can be generated by 
the Tayler-Spruit dynamo. 

of the convective core, the dependence of the coefficient 
-Dchcm, in the power 6 of {-^) reduces the magnetic dif- 
fusion drastically, so that as mentioned in Sect. 4.2 the ra- 
tio of the transport of chemical elements by the magnetic 
instability to the transport by circulation may amount to 
about 1 order of magnitude. Some tests indicate that this 
makes the chemical enrichments in helium and nitrogen 
at the stellar surface are stronger, but not too different, 
from those without magnetic field, despite the fact that 
the diffusion coefficients with magnetic field are orders of 
magnitude larger over most of the stellar interior. On the 
whole, we see that for chemical mixing also, the magnetic 
instability plays a great role. 

It is also interesting to see that diffusion coefficients by 
the magnetic field are almost equal (transport of chemical 
elements) or even larger (transport of the angular mo- 
mentum) than the thermal diffusivity. This means that 
magnetic effects by the Tayler-Spruit dynamo are in gen- 
eral equal or larger than thermal effects. Globally thermal 
effects may be relatively more significant in the outer lay- 
ers. 

6.2. Existence of a magnetic field due to tlie energy 
condition of differential rotation 

An important question in the models is to determine at 
each shell mass whether differential rotation described 
by parameter q is sufficient to create the Tayler-Spruit 
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Table 1. Structural parameters of the model of 15 M0 with Vi^i = 300 km/s when = 0.571 



M./Mq r/Ro n Vj, V,ad Vad g S Hp K 

6.5726 1.445 0.76E-04 -O.llE+00 0.18E+00 0.317 0.337 0.86E+05 1.334 0.33E+11 0.35E+10 

11.0386 2.080 0.69E-04 -0.33E+00 0.32E-03 0.269 0.348 0.69E+05 1.255 0.27E+11 0.79E+10 

14.9421 4.159 0.60E-04 -0.84E-01 0.14E-06 0.225 0.326 0.23E+05 1.455 0.19E+11 0.81E+12 



Table 2. Diffusion coefficients of the model of 15 Mq with Vmi = 300 km/s when Xc = 0.571 



6.5726 
11.0386 
14.9421 



0.48E-06 



0.73E-07 



0.84E-09 0.26E-06 



0.17E-12 0.17E-06 



0.12E+13 
0.51E+19 
0.16E+26 



0.17E+09 
-0.66E+01 

0.31E+19 
-0.57E+08 

0.24E+28 
-0.23E+14 



i^angl 

0.15E+14 
0.14E+14 



0.68E+10 
-0.77E+02 

0.15E+11 
-0.16E+03 



0.29E+15 0.19E+12 



.DanglP 

0.55E+14 
0.52E+14 



0.15E+14 
-0.42E+03 



0.49E+11 
-0.29E+03 

O.lOE+12 
-0.59E+03 

0.23E+10 
-0.23E+02 



IJshcar 

0.40E+04 
0.28E+07 
0.19E+09 



C/circ 

-0.19E-03 
-0.23E-04 
0.13E-03 



0.1 



0.05 




Fig. 8. The quantity - y-j^ j y-p-j^^ ^ 
Eq. (|39|l vs. r/RQ. When this quantity is positive in the 
region where case 1 applies, i.e. above r/R© = 2.69, differ- 
ential rotation is sufficient to generate the magnetic field, 
which is the case here. 

dynamo to produce a magnetic field. We examine here 
whether these conditions are fulfilled in the model with 
rotation only studied in the previous subsection. Firstly, 
we examine the regions adjacent to the core between 
r/R0 = 1.32 and 1.59, where case applies, since iV^ dom- 



(JL\ 



( ^ V 



given by 



inates. Fig. [3 shows the difference we see that 

in the concerned region between t/Rq = 1.32 and 1.59, 
this difference is negative, thus magnetic field is present 
there. This is an interesting result because it means that 
despite the strong restoring buoyancy force due to the very 
large ^-gradient, the differential rotation is high enough 
to develop magnetic instability. At each time step during 
evolution, such tests need to be performed. 

Secondly, we examine the zone above r/R© — 1.59 and 
up to r/R© = 2.69, where case IP applies. There, we check 
for the difference — If this expression is negative, 
magnetic field is created. This expression lies between - 
0.40 and -0.50 in this whole intermediate region. Thus, we 
conclude that magnetic field is also present there. 

The criterion for the existence of the magnetic field in 
the external zone, which corresponds to case 1 is given by 
Eq. H39|l . From Fig. ©, we see that in this zone which 
lies between r/R© — 2.70 and the surface, the differ- 
ence ^ > * ( r'^N-T ) * generally positive, which 
means that magnetic field is present over that region. 
However, we notice that very close to the surface this dif- 
ference goes to zero. This means that at the surface, dif- 
ferential rotation is becoming insufficient to generate the 
magnetic instability. This is interesting because it may 
explain why there is in general no strong magnetic field 
observed at the surface of OB stars (cf. Mathvs l^()()3|l . de- 
spite the likely existence of a strong internal field created 
the Tayler-Spruit instability. 

These results show that in a rotating star, the con- 
ditions on the differential rotation for the growth of the 
magnetic field are largely realized, with the possible ex- 
ception of the superficial layers. 
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7. Conclusions 

The main conclusion is that the Tayler instability and the 
Tayler-Spruit dynamo are of major importance for stellar 
evolution, both for the transport of angular momentum 
and for the transport of chemical elements. Future evolu- 
tionary models applying the results of this work will be 
made to study the results on tracks, surface composition, 
rotation, etc... 

It is likely that in a rotating star calculated with mag- 
netic field from the beginning, the differential rotation is 
very much reduced by the magnetic transport of angular 
momentum described above. We may suspect that differ- 
ential rotation will be reduced down to a stage where the 
criterion H37|l discussed in Sect. 5.1 is just marginally sat- 
isfied, i.e 

1^1 - ^ • (47) 

Indeed, if differential rotation is higher than given by this 
criterion, magnetic field develops and the associated cou- 
pling reduces differential rotation. If, at the opposite, dif- 
ferential rotation is lower than given by criterion H47|l . 
the growth of mean molecular weight by nuclear reactions 
in the central regions together with angular momentum 
conservation will produce an enhancement of differential 
rotation. Thus, a stage of marginal equilibrium is most 
likely reached during MS evolution. It may also be that 
the outer layers never have sufficient differential rotation 
to build Tayler-Spruit dynamo. Further numerical models 
will explore the evolution of stars with rotation and mag- 
netic field, and analyse the coupling between magnetic 
field and differential rotation. 
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